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Abstract 



We study the statistical properties of the spectrum of a quantum dynam- 
ical system whose classical counterpart has a mixed phase space structure 
consisting of two regular regions separated by a chaotical one. We make 
use of a simple symmetry of the system to separate the eigenstates of the 
time-evolution operator into two classes in agreement with the Percival clas- 
sification scheme ||. We then use a method firstly developed by Bohigas et. 
al. || to evaluate the fractional measure of states belonging to the regular 
class, and finally present the level spacings statistics for each class. The level 
spacings distribution of states belonging to the irregular part of the spectra 
as well as that of the complete set of levels corroborate the Berry-Robnik 
surmise ||. We further present a statistical study of the regular levels. The 
presence of intermediate states - states which belong to neither class as long 
as h is finite, phase spatially mixed among the set of regular ones, together 
with the small fractional measure of regular states strongly affects the corre- 
sponding level spacings statistics, resulting in a non negligible deviation from 
the expected Poisson distribution. We see however the remarkable agreement 
of the irregular level spacings statistics as a direct confirmation of the Berry- 
Robnik surmised. 
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1 Introduction 



For more than a decade, the study of quantum mechanical systems whose 
classical counterpart exhibits chaos has attracted much interest. One moti- 
vation for this study is the paradoxical fact that while the correspondence 
principle, as we understand it, should imply a quantum manifestation of clas- 
sical chaos, the Schrodinger equation is linear. As a consequence, the time- 
evolution operator is unitary, and this suppresses any exponential divergence 
in the time evolution of quantum states. As a spectacular manifestation of 
this fact, time-reversal invariant models show no lost of memory : reversing 
the time at a certain moment T brings us back to the initial situation after 
another time interval T, while this would require infinite precision in a clas- 
sical chaotic system. Thus a basic manifestation of classical chaos seems to 
have no place in quantum mechanics. 

On the other hand, the destruction of an integral of motion, of a quantum 
number, has striking effects on the statistical properties of quantum spectras. 
It is today taken as granted that in a classically integrable system, the levels 
are uncorrelated, and so have a poissonian level spacings distribution || (a re- 
markable exception being the one-dimensional harmonic oscillator) and that 
in classically fully chaotic models, the level spacings distribution has a dra- 
matically different shape : it obeys predictions of random matrix theory, i.e. 
it exhibits level repulsion HJ. The situation in mixed systems, where regular 
and chaotic regions coexist in the classical phase space, is more intricated. 
In an old paper Percival [[| classified the eigenf unctions of the Schrodinger 
equation into two classes belonging to either the regular regions, where the 
invariant tori are not destroyed, or the chaotic one. This classification was 
based mostly on the correspondence principle and has been numerically con- 
firmed a few years ago by Bohigas et. al. 0. While the eigenf unctions 
that are mostly confined on classically regular regions - we will call them the 
regular eigenf unctions - tend to concentrate on invariant tori, the irregular 
ones tend to spread uniformly over the chaotic region as h — > 0, as has been 
rigorously demonstrated by Shnirelman This picture is assumed to re- 
flect reality in the semiclassical limit h = 0. Following this classification, 
Berry and Robnik postulated that the part of the spectrum that corresponds 
to regular eigenfunctions, has a poissonian level spacings distribution in op- 
position to the one corresponding to the irregular eigenstates which exhibits 
level repulsion || . This surmise led them to an expression for the level spac- 



3 



ings distribution for mixed systems that has been observed convincingly only 
recently for the case of the kicked rotator on a torus. As pointed out by 
Prosen & Robnik and Li & Robnik || |TB|, reasons for this difficulty of ob- 



servation could be that we are not deep enough in the semi-classical regime. 
As long as h is finite, a certain number of wave-functions belong neither to 
the regular nor to the irregular set of eigenfunctions. We may think of states 
making use of the Heisenberg uncertainty to overlap the frontier between the 
regular and irregular regions of the classical phase space, or states located 
on the regular region which, due to the finiteness of the Planck constant, do 
not yet belong to the set of regular states. Consequently, the Berry-Robnik 
regime should be observable only in the far semi classical limit. We will come 
back to this point later. 

In this paper we present a spin model allowing a precise study of a mixed 
regime. The reasons for this are first that, in a special regime, an approx- 
imate simple symmetry of the phase space structure, namely S z — > —S z , 
allows the separation of regular states from the irregular ones, and secondly 
that the frontier between the regular and the chaotic zones is rather sharp, 
thus minimizing the number of intermediate eigenstates. This enables us to 
compute the level spacings statistics independently for the regular and irreg- 
ular states. We interpret the fact that these statistics obey quite well the 
poissonian distribution and the GOE respectively as a direct confirmation of 
the Berry-Robnik surmise. 

We study the quantum system defined by the following Hamiltonian |: 

U gm : = |((1 - z 2 )S 2 z - z 2 Sl) + kS z A t (1) 
and the corresponding unitary time evolution operator : 

Ut := exp(-i«S,) «p(~y((l - z 2 )S 2 z - z^ x )T) (2) 

where S = (S x , S y , S z ) = h(s x , s y , s z ) = hs are spin operators satisfying the 
usual commutation rules (e 4Jlfc is the total antisymmetric tensor of 3 rd order): 

[S^S,] =ihe^ k S k (3) 



4 We use bold characters in the quantum case in contrast to normal ones which refer 
to the classical and semiclassical cases : S refers to the quantum spin operator while S is 
either a classical or a semiclassical spin. 



4 



< k < 2tt, A t := En=-oo 5 (t ~ nT ), N = energy" 1 time" 2 and < z < 1 
. Models of this kind have been extensively studied |ll|. They represent a 



spin which evolves under the influence of a classically integrable Hamiltonian 
H° m = y ((1 — z 2 )S 2 z — z 2 S 2 x ) during a time T after which the spin undergoes 
a rotation of angle k around the z-axis. The regime we consider is defined by 
k = 1.1, T — ^ and z 2 = |. Classically there are two regular zones around 
the north and south poles surrounding a chaotic region which is fairly well 
symmetric under S 2 reflection (Fig.l). In the semiclassical limit which cor- 
responds to hs = S = constant, h = s^ 1 — > 0, states which are located on 
the chaotic region tend to cover it homogeneously according to Shnirelman's 



theorem [12j]. Since this region is symmetric under S z reflection, the expec- 
tation value < ^chaos\ s z\^ chaos > of such a state tends to disappear as we 
approach the semi-classical limit. For small but finite h, the distribution of 
< ^fcls^l^/fc >, where > is an eigenstate of the operator Ut defined in 
(2), will then present a sharp peak around zero corresponding to the irregu- 
lar states surrounded by two smaller bumps corresponding to regular states 
(Fig. 2). This allows us to separate easily the regular states from the irregular 
ones, the validity of this selection being confirmed by a numerical semiclas- 
sical argument presented in section 3 as well as an extensive study of the 
Husimi distributions of the selected states [l^ . 



The paper is organized as follows : Section 2 is devoted to a short presenta- 
tion of the classical model. In section 3 we derive some useful semiclassical 
quantities such as the density of states and the expression for the action. 
This will allow us to estimate the number of regular states, and give a check 
of our selection criterion. In section 4 we present the quantum mechanical 
model as well as our numerical results for a spin magnitude s=500. All of 
them were obtained using direct diagonalization techniques. Conclusions and 
further remarks are given in section 5. 



2 Classical model 

The unperturbed classical Hamiltonian 

AS := " ^ " ^) (4) 

has two degrees of freedom and is an integral of motion. The trajectories are 
confined to the intersections of the sphere |o| = S with the cones of constant 
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energy E = f ((1 - z 2 )S 2 z - z 2 S 2 x ). The perturbation 

Hh := kS z A t 



(5) 



corresponds to a rotation of angle k around the z-axis performed at time 
intervals T. Its addition leads to the destruction of the energy surfaces, 
and allows more and more trajectories to wander chaotically on the sphere of 
constant spin magnitude as k and T grow. Expanding up to the first order 
in SS Z :— S — S z near the poles S z = ±S, we get a one-dimensional harmonic 
oscillator of period T = 2\[2^. In particular we have 8S Z = 0(5S 2 ) : in this 
approximation 5S Z is an integral of motion and is furthermore conserved by 
the perturbation H\ x too. It is thus conceivable that the invariant torii near 
the poles will offer more resistance to the perturbation than those located 
away from them. We use this property to find a regime in which there are 
two regular islands around the poles approximately related by the operation 
S z — > —S z and separated by a chaotic region. This we achieved by setting 
k = 1.1, T = -^g,z 2 = 0.5 (Fig.l). The regular islands occupy in a good 
approximation the region 0.22S' 2 < E < E max = 0.25S' 2 . 



3 Semiclassical approach 

We compute the Green function for a trajectory of positive energy and the 
density of states for the unperturbed case T = z 2 = 0.5. We follow the 
lines drawn in |§. We first write the unperturbed Hamiltonian in canonical 
variables (S z , <fi) for the chosen regime : 



m 



H = ^(^(l + cos 2 (0))-5 2 cos 2 (0)) 



(6) 



The action integral for a trajectory of energy E starting at 0o and ending at 



reads (S = hs, e 



h m ' 



s, 



00 



4e + s 2 cos 2 ( 



cos z 



2tt 



4e + s 2 cos 2 ( 



o \ 1 + cos 2 0') 



(7) 



where we have set <fi* = <fio + ((f) — <po) mod 2-k, and np is the number of 
complete revolutions accomplished between O an d <p (4> = 2imp + (f>*). The 
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sum runs over classical orbits (3 of constant energy. This leads us to the 
expression for the corresponding Green function : 



g(<f>o,(j>,e) = --^ v / |detPi i/3 ((/)o,0)|exp(- l S / 3((/)o,0,e) - — 



with 



det£> li/3 Oo,0) = 



g)W 
(E — E n + ie) 

d 2 S p d 2 Sf3 d 2 S p d 2 Sf3 



(8) 



'o de 2 d(J)de d(ft de / 
-4 



m 2 y / (4e + s 2 cos 2 (0 o ))(l + cos 2 (0 o ))(4e + s 2 cos 2 (0))(l + cos 2 (0)) 



(9) 



this result being obtained by partial differentiations of (7). Since this latter 
value never changes sign, the Maslov index lp vanishes and so the diver- 
gence of the Green function leads to the following semiclassical quantization 
condition : 



2n 



4e + s 2 cos 2 ( 



o \ 1 + cos 2 (0') 



-d(j>' = 2ttM 



(10) 



for any integer < M < s. We have then for the averaged density of states 
(N is the number of states) : 



P(e) = 



1 



-Im 

Ntt 

1 

irh Jo 

2 rl-K 



2tt 



d(f>g(<j),<j),e) 



#^|det 

d(j) 

Jo 



Jo ^{Ae + s 2 cos 2 (</>))(1 + cos 2 (</>))) 



(11) 



The last equation states in particular that the averaged density of states is 
proportional to the classical orbit period. Fig. 3 shows the agreement of this 
semiclassical result with the numerically obtained density of states for the 
unperturbed quantum model at s=1000. Using (11) we estimate the number 
of states occupying the regular region of figure 1 : 

r0.25s 2 



N reg ~ / p(e)de 



1 

24 



(2s + 1) 



(12) 
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The number of states occupying this region in absence of perturbation is ^ 
times the total number of states. This gives us a first approximation for the 
number of regular states we must select. A better approximation in presence 
of perturbation is given using a method developed by Bohigas et. al. 0. We 
must evaluate the number N reg of trajectories that satisfy the condition : 



E 

i=0 



2ttM (13) 



for some integers M and P, while it nearly closes on itself after the P th kick, 
i.e. : s^Xj) ~ s z ((f>Q ). <p~ is the angle between the x and the y component 
of the spin just before the i th kick while <pf refers to the same angle right 
after this kick. 

This condition means that the action integral must still be an integer multiple 
of 2 7r, and that simultaneously, the orbit must be closed. This condition has 
meaning only on regular regions were the invariant torii are not destroyed, 
so that the integrals make sense. We transform this condition and compute 
the number of trajectories satisfying : 



fl +1 S z {(j)')d(j)' + S z ((j)j)K 



M (14) 



for integers M and P, and P sufficiently large. With this we replace two 
conditions by only one numerically more tractable condition. Since our task 
is to evaluate the number of regular semiclassical levels, and not to determine 
them precisely, we believe that condition (14) is sufficient. The number of 
regular states we numerically estimated with (14) is 50 ± 4 for s=500, i.e. 
slightly larger than that estimated with (12). In the next section, we will 
consider this estimated number of regular states as a check of the validity of 
our selection criterion. 



4 Quantum Model 

In this section we study the statistical properties of the spectrum of the 
quantum Hamiltonian (1) for integer spin magnitude. Since the perturbation 
term is time-dependent, the energy is no longer a good quantum number, 



S 



and we are led to define quasi-energies and quasi-energy eigenstates. The 
Schrodinger equation leads to the following time evolution from right after a 
kick to right after the next one: 

vf(T+) = U T M>(0 + ) = exp(-i«S,) exp(~H° qm T)^(0 + ) (15) 

Quasi-energies A and quasi-energy eigenstates ^\ are then defined by : 

U r tf A = expHA)^ A (16) 

Since Ut is unitary, the A's are real and defined modulo 2 it. We introduce 
two parities : 

U\fi >= | - fi > (17) 

e\fi >= (-iy-^\n > (is) 

We can express the time-reversal operator T in term of these two parity 
operators : 

no @\n >= T\/j >= (-1) S ~ M | -fi > (19) 

In the integer spin case the eigenstates > of Uj. := exp(— |H° m T) satisfy 
the conditions : 

n|^>=±|#> (20) 
0|^>=±|^> (21) 

So H° m and are in particular time-reversible. The perturbation breaks 
the II-parity but leaves the 0-parity unbroken (]. We will concentrate on the 
study of even states, i.e. those states satisfying : 

0|tf >= |$ > (22) 

However partial results obtained for the odd set of states corroborate the 
results presented here. The key point is now to find a clear quantum mani- 
festation of the approximate symmetry S z — > — S z of the classical phase space 

5 In the half-integer spin case, eigenstates of the unperturbed time-evolution operator 
are eigenstates of the ©-parity only. The latter is left unbroken by the perturbation we 
consider. 
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structure. A practical solution is given by Shnirelman's theorem which states 
that in the semi classical limit, the quantum states that are confined on the 
classically chaotic region of the phase space tend to cover it uniformly. To 
get an insight in this statement we use the following resolution of unity fl4 



2s 



[ d9d<p sin 9\9,<l> >< 9,<p\ (23) 



7T 

where we introduced coherent states of the spin SU (2) group : 



\9,<P>:= ± /(^ysin(^)^cos(^ + ^-^|^> 

These are states that are centered on the point (9, (/)) of the sphere and which 
minimalize the quantum uncertainty. 9 is defined by S z = Scos(9). Using 



(22), the symmetry of the chaotical region and Shnirelman's theorem [12 



, ,„ . f on the regular region 

< ^l*. >^ { congt on the ch g aotic J gion (24) 

it is then easy to show that 

< y c haos\Sz\y chaos >^ (25) 

in the semiclassical limit. This translates into Fig. 2 where we plotted an his- 
togram of the expectation value of s z taken over quasi-energy eigenfunctions 
for s=500, k = 1.1, z 2 =0.5, and m=l. The central peak clearly reflects our 
reasoning, while the two smaller bumps surrounding it are mainly due to the 
regular states that are confined to the classical stability islands. The gap 
in-between is a consequence of the uniform distribution of irregular states. 
It is remarkable that this gap overlaps the classical frontier between regular 
and chaotic region. 

We used this property to part the irregular states from the regular ones and 
then study separately the statistical properties of the spectrums of each class 
of states. We believe this criterion is justified since the fluctuations 



As, = ^< S 2 > - < s z > 2 (26) 

of regular states is much smaller than the "Shnirelman gap" appearing in 
the histogram of Fig. 2 between the huge central peak and the smaller 
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bumps. As a consequence only very few regular levels will be selected with 
the set of irregular ones, while maybe more irregular will be counted with the 
regular ones. Moreover, the fact that the number of selected regular states 
is in complete agreement with the numerical semiclassical evaluation given 
by (14) confirms the relevance of this selection criterion. 
We now turn our attention to the study of the spectral properties of the 
time-evolution operator (16). Due to the ©-symmetry (18), belongs to 
the circular orthogonal ensemble, and not to the circular unitary ensemble 
as would be expected from the fact that the perturbation breaks the time- 
reversal symmetry. This situation is similar to the one encountered by Berry 
& Robnik in certain Aharonov-Bohm billiards ||, or by Delande & Gay in 
the Hydrogen atom in a magnetic field |TD[ where the system violates the 



time-reversal symmetry, but possesses an invariance under a combination of 
the time-reversal and another symmetry, in our case the II-symmetry. We 
thus expect a linear repulsion for the part of the spectrum belonging to the 
irregular states. 

The results of our study for a spin magnitude s=500 are plotted in Fig. 4 
to 9. Fig. 4 shows a plot of the level spacings statistics for 4233 irregular 
level spacings computed by diagonalizing ten different evolution matrices for 
T = and 1.05 < k < 1.15. The solid line is the predicted Wigner distri- 
bution. The agreement is excellent. In Fig. 5 we plotted the corresponding 
cumulative level spacings distribution defined in term of the level spacings 
distribution P(s) by 

W(s) = !~dtP{t) (27) 



o 

Also shown are the Poisson and the Wigner distributions. As shown in inset, 
small deviations from the Wigner distribution appear only around spacings 
s ~2, but have no significance to our opinion. 

Fig. 6 and 7 show the level spacings and the cumulative level spacings distri- 
bution for 572 regular levels taken from twenty different evolution matrices 
for T = -i^ and 1.05 < k < 1.15. The difficulty here is the relatively small 
number of regular states (~ 25). Accordingly, only few intermediate or irreg- 
ular states can have relatively big effects on the statistics. In a convenient 



D We have used s for the level spacings to avoid confusion with the spin magnitude. 
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basis, Ut can be represented by the following matrix : 



U T : 



( 



R K 
K C 



) 



(28) 



R = R{5i : j is a diagonal N reg x N re5 matrix which corresponds to the N re5 reg- 
ular states, C = CiSij is a diagonal N c haos x N chaos matrix which corresponds 
to the N c haos irregular states, and K = 0(h) couples the two subspaces as 
long as h is finite. In our picture K disappears in the semiclassical limit, 
and the Ri and satisfy a poissonian statistics and a GOE statistics re- 
spectively It would be of course a hopeless task to try to determine K for 
finite h. The important point is to recognize that as long as h is finite but 
small enough, K couples only few regular states with irregular ones, this fact 
resulting in a deviation from the Berry-Robnik surmise. This deviation is 
then naturally much more important for the regular part of the spectrum, 
since it contains much fewer levels than the irregular part. We believe that 
this is the reason for the deviation of the statistics of the set of levels we have 
selected as regular from the poissonian predicted behaviour. We must recall 
that our whole reasoning is based on the assumption of two classically homo- 
geneous stability islands. In such a case, semiclassical wave-function would 
mimic classical orbits and would therefore fit together as concentric circles. 
The presence of hyperbolic fixed points or cantori may change this picture, 
possibly turning regular states into intermediate ones as long as % is finite. 
The semiclassical wave-function overlap and thus interact at certain regions, 
and this, in the Pechukas picture modify very sensibly the equations of 
motions governing the evolution of the quasi-energies A as k or T is mod- 
ified, resulting in the appearance of level repulsion. So some intermediate 
states are phase spatially mixed among the set of states we have selected as 
regular and consequently modify the corresponding statistics. Their effect is 
furthermore enhanced by the small ratio of regular levels. A current investi- 
gation of the Husimi densities of the selected regular states corroborates this 
reasoning ||15|| . Finally we show in Fig.8 and 9 level spacings and cumulative 
level spacings statistics for the complete set of levels. We compare our results 
with the Berry-Robnik prediction for a fractional measure of regular states 
as approximated by (12). The agreement is amazing, and corroborates our 
picture. The x 2_ test for both graphs - x 2 =25, i.e. half the number of boxes 
for Fig.8, and x 2 =1480, i.e. 3.3 times less than the number of levels for Fig. 9 
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- gives full statistical significance to these last graphs. We see them as a good 
evidence for the validity of the Berry-Robnik surmise in our model. 

5 Conclusion 

We studied the statistical properties of a quantum spin model whose classi- 
cal counterpart exhibits a mixed phase space configuration. Due to a simple 
approximate symmetry, whose effect on the quantum system is drastically 
enhanced by Shnirelman's theorem, we were able to separate the irregular 
from the regular levels, thereby confirming implicitely the validity of the Per- 
cival classification. We then performed a separated statistical study of these 
levels. The results confirm the Berry-Robnik surmise : while the irregular 
set of quasienergies exhibits a clear wigner-like shape, the regular part of 
the spectrum has a clearly different shape, though its spacings distribution 
does not follow strictly a poissonian law. This deviation is interpreted as 
the presence of both irregular and intermediate states among the selected 
regular ones, their effect being enhanced by the relatively small number of 
the latter. Nevertheless, due to the small number of regular states we believe 
that the irregular statistics is much more significant, and see our results as a 
good confirmation of the validity of the Berry-Robnik surmise in our model. 
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Figure Captions 



Fig.l: Orthogonal projection of the classical phase space on the (S x ,S y ) 
plane for the case T=^|, k = 1.1 and z 2 =0.5. 

Fig. 2: Histogram of the expectation value of S z taken over the eigenstates 
of the unitary time evolution operator defined in (2). 

Fig. 3: Density of states for the unperturbed Hamiltonian according to (11) 
(solid line) as compared to numerically obtained datas for the case s=1000 
(squares). 

Fig. 4: Level spacings distribution for 4233 irregular level spacings obtained 
through direct diagonalization of ten evolution matrices in the parameter 
range T = ^ and 1.05 < k < 1.15. 

Fig. 5: Cumulative level spacings distribution for the same case as fig. 4. 
In inset : regions of small deviation relatively to the Wigner-distribution. 

Fig. 6: Level spacings distribution for a set of 472 regular level spacings 
obtained through direct diagonalization of twenty evolution matrices in the 
parameter range T = and 1.095 < k < 1.105. The solid line is the pre- 
dicted Poisson distribution. 

Fig. 7: Cumulative level spacings distribution for the same levels as Fig. 6 
compared to the Poisson distribution. 

Fig. 8: Level spacings distribution for a set of 5000 regular and irregular 
level spacings obtained through direct diagonalization of ten evolution ma- 
trices in the parameter range T = ^ and 1.095 < k < 1.105. The solid line 
is the predicted Berry-Robnik distribution with fractional measure of regular 
states p\ = 0.08. x 2= 25 is half the number of boxes. 

Fig. 9: Cumulative level spacings distribution for the same levels as Fig. 8 
compared to the poissonian and the Berry-Robnik predicted distribution. In 
inset : Same curve compared to the Wigner distribution. x 2 =1480 is 3.3 
times less than the number of levels. 
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